Softening and melting of a vortex lattice in presence of point disorder. 



OO 



Denis Feinberg 

Lahoratoire d'Etudes des Proprietes Electroniques des Sohdes, Centre National de la Reeherche Scientifique, 

BP 166 38042 Grenoble Cedex France 
(February 1, 2008) 

A phenomenological model is proposed for melting of a vortex lattice, based on screening of 
the elastic shear modulus by mobile or partially pinned dislocations. A first-order softening line 
is found and ends at a critical point beyond which the lattice crosses over to an hexatic vortex 
solid. The consequences of softening on vortex dynamics are explored, as fingerprints of plastic 
dynamics: a reentrance of single vortex behaviour, for both depinning and collective creep, occurs 
as the field increases, with non-monotonous creep exponents. This general scenario is supported by 
recent experiments in high- Tc materials and suggests that for a 3D vortex lattice at low temperature 
the field induces a continuous order-disorder transition towards a glassy phase. 
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PACS numbers: 74.60. Ge, 64.70.Dv, 62.20.Fe 



The interplay between thermal fluctuations and 
quenched material disorder controls the phase diagram of 
vortex matter in High Temperature Superconductors [Q . 
In clean crystals, a first-order melting line exists but ends 
at a critical point , while in strong fields or for strong 
enough material disorder, the vortex lattice (VL) melts 
continuously, possibly through a vortex glass (VG) tran- 
sition 0] . It has been recently proposed that a soHd-solid 
transition also occurs with increasing field, between a 
quasi-ordered "Bragg glass" and a disordered vortex solid 
[^-0. This is supported by irreversible magnetization 
measurements showing an increase of the effective criti- 
cal current ( "second peak" ) and by neutron diffrac- 
tion experiments [p^ . Such a transition is attributed to 
the proliferation of disorder-induced topological defects. 
However, contrarily to the pure melting transition, its 
first-order character might be questioned, since those de- 
fects are also pinned by point disorder and therefore will 
not dramatically screen the elastic shear modulus, con- 
trarily to a pure melting scenario. In what follows we 
propose a simple model of plastic behaviour, based on 
short range fluctuations. It produces a first-order transi- 
tion with a critical point, by capturing the two effects of 
point disorder, e.g. promoting and pinning dislocations. 
From the calculated renormalized shear modulus, vortex 
dynamics can be discussed. 

Let us consider the softening of the shear modulus 
Gee = (s^xy^ i4>o = ^ is the London length with field 



parallel to the c-axis for HTSC). For simplicity the T- 
dependence of Cee is neglected. As shown by Brandt 
and Marchetti and Nelson softening is due to 
motion of edge dislocation (DL) segments normally to 
the field. These are part of DL loops with screw com- 
ponents normal to the field. Their natural unit length 
is Ln ^ (■FT-Y^'^ the elastic intervortex "interac- 

tion length" ||l|]. Here cq and 7 are the vortex spacing 
and the anisotropy coefficient respectively and eo is the 
vortex line stiffness. Rather than a detailed description 



of these loops, let us forget their line character and at- 
tempt at a mean-field picture, valid in principle for high 
DL densities . First, the edge and screw core energies 
being roughly equivalent, on can simply calculate the en- 
ergy and average density of screw DL cores. Those can 
be considered as "kinks" where a vortex line passes from 
a valley to another of a periodic "egg-carton" potential 
representing the reference perfect VL on top of which 
DL's are created. This extended potential should replace 
the parabolic "cage" potential, valid only for the elas- 
tic analysis of the fluctuations of one vortex line in the 
mean field of the others [^j^. In what follows, we shall 
more simply obtain the kink density from a line's fluctu- 
ations in a renormalized cage potential. The kink energy 
is Ek ^ otEei where Eei ^ Acj^CeeaQLo is the typical 
elastic cohesive energy (a > 1 is a numerical constant 
and Ci ^ 0.15 the Lindemann constant). On the kink 
length ^ Lq, a vortex line experiences gaussian lateral 



fluctuations < 



^^(T + A) where A - Uq{^J^<- 
is the elastic pinning energy on the length Lq Here 
Lc and Uo are the line's Larkin length and pinning energy 
respectively. We shall consider here only independently 
pinned vortices, i.e. the single vortex (SV) regime for 
which Lq> Lc, then C, ~ 3/5 0]. This allows to write the 



DL density as Ud 



exp(- 



aE^i 
"T+A 



where Li 



is a longitudinal "domain size", and Ri ^ ao exp( j^_^_^ ) 
is a transverse positional correlation length. The picture 
here is that of an hexatic VL since DL's preserve orien- 
tational order. 

Let us now turn to the screening phenomenon. If edge 
DL's were free to move, they would perfectly relax shear 
stresses and renormalize Cge to zero p^ . However, they 
experience two types of potentials. The first one is due 
to the VL periodicity and builds the Peierls-Nabarro bar- 
riers, of order vE^i {v of order 1). The second one is the 
disorder pinning potential, which in the SV regime acts 
on DL cores of length Lq , thus is also of order A and has 
a correlation length Sq ~ cq. In the bundle or 3D regime 
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{Lc > Lq) the situation is more complicated : the part of 
the DL pinning energy coming from elastic displacements 
of neighbouring vortices builds a random potential which 
has a priori a correlation length Sq > oq. For a dense 
collection of DL's, DL interactions will lead to collective 
pinning effects, e.g. new pinning lengths will appear. 
Such" plastic" correlation lengths should be considered 
in a refined analysis. 

In the above picture, the energy and length scales are 
the same as in the elastic analysis. Nevertheless, the 
self-consistent screening mechanism leads to qualitatively 
different results for the melting line. To analyse the soft- 
ening of the vortex solid by proliferating DL's, let us 
define a renormalized shear modulus C|g = xCee{x < 1). 
Then all relevant "elastic" quantities must be accordingly 
modified as L§ = x-^/'^Lq > Lq, = E^ix^/"^ < E^i 
and — Aa;^^/^° > A. Softening therefore increases 
the density of DL's, which in absence of disorder drives 
the first-order melting transition. To close the self- 
consistency loop one must express the softening of Cge 
by moving DL's. This can be done similarly to Ref. 
though in a simplified way. Indeed, vortex line displace- 
ments are planar and one can write a simplified screening 
equation for each "layer" of vortex matter of thickness 
L^, e.g. ^ = ^ X where x = ^NdL^' < >d is a 

"plastic" susceptibility due to motion of Nd ~ exp(^^) 
edge DL segments of length L^. The correlation func- 
tion < it^ > of DL displacements can be evaluated, 
given the potential well (Peierls-Nabarro plus disorder) 
in which they sit. One can write < >~ cioUfiv) with 
y — The function f{y) describes the depin- 

ning phenomenon : equal to one when y << 1 (strongly 
pinned DL), it increases rapidly when y > 1. A simple 
choice is here f{y) ~ exp(y). Other choices lead to qual- 
itatively similar results. 

It is essential to include the effect of thermal depin- 
ning. This can be easily done within the SV regime 
by taking L,(T) ^ L^(o)^ exp((^)3) and A(r) ^ 

A(0)(^|^)-6/5 if T > C/o §. This, together with the 
above definitions and with the use of renormalized quan- 
tities, leads to the self-consistency equation 



A similar expression could be written in the bundle 
regime, in which one notices that the plastic suscepti- 
bility is larger by a factor (f^)^ due to weaker DL pin- 
ning. The solution of Eq. (1) is obtained as a function of 
T = and S = ^ ' which control the importance of 
thermal and disorder fluctuations respectively. First, in 
the case of zero disorder, one sees from Eq.(l) that the 
screening effect becomes dramatic when T ~ E^^, which 
determines the "pure" melting temperature. Here the 
constants u and a are chosen so as to yield the melting 



transition at = E^i- We take v = 0.5 and a = 8. The 
transverse correlation length Ri{t) is plotted on Fig.l. 
At the transition this corresponds to a correlation length 
of order 10^ unit cells e.g. a quasi-perfect lattice {x ~ 0.8, 
a value depending on the choice of a). 

For non-zero disorder we take as a reference the field 
Be{0) which marks at T = the transition to the disor- 
dered phase in the Lindemann-based analysis . Then 
one defines b — ^ ^ 6^^"^ (SV regime) , and similarly 

Te = T/Te = where = T^{B,) (see Fig.3). 

First, for weak disorder, the transition remains first or- 
der but occurs at Tm < T^, in qualitative agreement 
with experiments (|J|l and theory |jl||. At T the 
correlation length Ri is still very large. On the other 
hand, for strong enough disorder (S > 1.7), the tran- 
sition becomes a smooth crossover. This allows to de- 
fine a critical point {Bcr,Tcr) ending the first-order line 
(Fig.3). Taking for instance Uq = 0.8Te one finds that 
Tcr ~ 0.52Te and B^r ^ IMB^ , with x^r 0.28. This 
corresponds to a rather small Ri ~ 5ao. The transition 
can also be represented by plotting Cqq ~ Bx{B) as a 
function of i? at different temperatures (Fig.3). It is lin- 
ear at low fields and low temperatures, then bends down 
and sharply drops at the melting field B„i{T) if T > Tcr 
(first-order transition) or goes smoothly through a max- 
imum if r < Tcr- Physically, for T < Tcr DL's are more 
strongly pinned thus not mobile enough to drive a sharp 
collapse of the shear modulus. However, even limited 
DL fluctuations (or order oq) can have a sizeable screen- 
ing effect. Provided the drop of Cee is identified with a 
melting line, our model is in qualitative agreement with 
experiments |^,^. The case of Bi2Sr2CaCu20% is spe- 
cial, since experiments provide stronger evidence for a 
first-order field-induced transition at low T . However it 
might also involve a dimensionnal crossover, which is not 
present in our 3D model. 

Let us now analyse the consequences of the present sce- 
nario on vortex dynamics. First of all, it means that DL's 
also move to screen elastic stresses, e.g. the dynamics is 
in reality plastic. A refined theory should take a detailed 
account of the dynamics of both elastic domains and DL's 
as a function of the length scale L{J) involved at a given 
current J. Here we simply analyze critical currents and 
collective creep (CC) within the renormalized elastic pic- 
ture, and focus on the boundary between SV and 3D be- 
haviours, defined by L{J) ~ Lq. It depends crucially 
on the vortex interactions thus on the effective shear 
modulus. Here, softening clearly favours a SV regime 
by decreasing the vortex correlations. First, for the SV 
critical current J^^,, L{J^^) = Lc and 3D pinning occurs 
for Bx{B) > Bsi, ~ ^^^2 ■ Since Bx{B) has a maxi- 
mum (see Fig.3), two situations are possible, depending 
on whether max(J5a;(-B)) is larger or smaller than B^b- In 
the first case, occurring at low temperature, a 3D regime 
occurs followed by a reentrant SV regime at high fields 
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(Fig. 2). In this case x{B,T) should in principle be re- 
calculated in the 3D regime by generalizing Eq.(l), but 
this would not change the boundaries of the SV regime 
(Fig. 3). In the second case, at higher temperatures, the 
bundle pinning regime is never reached, vortices remain 
pinned individually. Then one must look for a bundle col- 
lective creep (CC) regime, and examine the frontier be- 
tween SV and bundle CC. It occurs at a current density 
Jsb{B) - Jsvi^y^^" 0- As before, a reentrant bound- 
ary is obtained by using the "rescaled" field Bx{B, T) 
instead of B, and is plotted on Figs. 2 and 3. While a reen- 
trance of SV critical currents was qualitatively suggested 
in the past as an explanation for the "peak effect" |^,^ , 
we here underline a similar phenomenon for creep, and 
provide a quantitative model. It is also much general 
since for small enough currents, in absence of screening, 
the dynamics would always be 3D. This has important 
consequences on CC behaviour. Indeed, in the 3D CC 
regime, if elastic^ vortex correlations result in a faster di- 
vergence of the creep barriers as J goes to zero : while 
the exponent ^ in U{J) ^ Uo{^Y is small, fisv ^ \ 
the SV regime, it is larger, | < < | in the elastic 
3D regime . This "elastic " analysis strongly disagrees 
with experiments performed in crystals of Y Ba2Cu^07 
0, [Nd, Ce)CuOi |l| and {Ba, K)BiOs. |6|jl|] as weU 
as HTSC films Q or multilayers Crystals show a 
"second peak" , together with non-monotonous creep ex- 
ponents, all features incompatible with the elastic analy- 
sis. Due to creep, the peak cannot be strictly attributed 
to a non- monotonous critical current, and a crossover in 
creep regime is possible However, this last scenario 
assumes bundle CC above the peak, with large /j,;,, while 
it was shown p6 17| that the CC scaling still holds pro- 



vided one uses a field-dependent ^i{B) which decreases 
at high fields. On the other hand, the reentrance found 
in our model nicely explains these results : as the field 
increases, the exponent /i has a maximum inside the bun- 
dle CC region and eventually goes back to the SV value 
as the SV regime is recovered. We thus claim that such a 
non-monotonous /i(i?) is a characteristic feature of "plas- 
tic" collective creep. Remarkably, in [Ba, K)BiO^ 
the dynamics above the second peak cross over to a VG 
transition, the dynamical exponent reaching a value close 
to the SV exponent /i^-u = y. Such small values seem 
ubiquitous close to VG transitions, an appear as a signa- 
ture of a renormalized SV creep regime. We thus show 
that softening is an essential phenomenon to account for 
the dynamics at and above the peak. More generally, 
the present scenario could also be applied to the criti- 
cal current "peak effect" pl|-p^ occurring close to i?c2, 
provided the full B— dependence of Cge is used. Recent 
neutron scattering experiments in Nb indeed show the 
succession of vortex lattice, hexatic and amorphous solid 
as B increases towards Bc2 |Q- 

The physical picture in the reentrant SV regime is the 



following : each vortex line wanders laterally on distances 
much larger than oq, thanks to edge DL's. Softening 
make the effective correlation between vortex lines much 
smaller, providing a simple mean-field picture of an en- 
tangled hexatic solid. However, it does not allow to fully 
describe the highly disordered phase. First, above the 
first-order drop (T > Tcr), Cjg is not exactly zero there- 
fore the phase from this point of view is still a solid. 
Similarly, at low temperatures and high field the shear 
modulus does not vanish, as expected close to a VG tran- 
sition. In fact, other topological defects such as vacan- 
cies or intersticials should be considered |2^j2^, also a 
more powerful approach is necessary to describe a tru- 
ely amorphous glassy phase with possibly second-order 
(solid-solid and solid-liquid) transitions, joining at the 
critical point (Bcr, ?cr)- Our approach, limited to screen- 
ing induced by short-range fiuctuations, should be im- 
proved to describe such transitions. One can neverthe- 
less try to extend the present screening idea to the VG 
transition proposed by Fisher First, in an elastic pic- 
ture, thermal depinning is characterized by the equality 
between T and the pinning energy gained at the scale 
of the temperature-dependent Larkin correlation length 
Lc{T). This length becomes very large at high temper- 
ature where disorder is effectively washed out. On the 
other hand, in the vortex glass transition picture, a glass 
correlation length (T) is defined Q . The behaviour at 
length scales L < is critical, e.g. characterized by the 
equality between T and the vortex glass energy fiuctua- 
tions. In the plastic vortex lattice picture, these fiuctua- 
tions involve the pinning energy of DL's. When the DL 
density is so large that the vortex solid becomes amor- 
phous, due to fiux cutting and other defects, a relevant 
length scale Cc{T) must be defined. In this phase the 
effective disorder should be renormalized down by large 
scale fluctuations. We suggest that through this phase Cc 
go smoothly towards the vortex glass correlation length 
^,g{T), with ^,g{T) ~ ^o{T,g - T)"'' [|. Thc abovc 
screening model can be simply modified by replacing 
by ^„g(r). Being critical, at the length T is of the 
order of the pinning energy therefore < >~ in the 
expression for the plastic susceptibility x- This, together 
with E^' ~ ~ T leads to - 1 t4-( — )^ which 

close to the VG transition amounts to ~ {Tyg — T)^''. 
Full softening is completed at the vortex glass transition 
where C§^ becomes zero and the vortex glass melts. 

In conclusion, we have proposed a model for the soft- 
ening and melting of a vortex lattice in presence of point 
disorder. We find a critical point ending the first-order 
line, and a continuous crossover towards an hexatic glassy 
phase at low temperature, associated with a maximum 
in the shear modulus. Softening of the vortex solid mani- 
fests itself in vortex decorrelation thus a reentrance of sin- 
gle vortex pinning and creep, providing an explanation 
for thc anomalous creep exponents. This model allows 
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to test quantitatively the possibility of disorder-induced 
plasticity, however much more work is required to analyze 
in detail the plastic dynamics and to properly describe 
the (amorphous) vortex glass transitions. 
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and increasing disorder, from top to bottom : 5 = 0,1, 1.71 



FIG. 2. Schematic {B, J) phase diagram showing the SV and 3D creep regime boundaries with renormalized Cee (thick lines) 
and the reentrance at high fields. Bi and B2 are the limits of SV critical currents. The dotted lines stand for the pure elastic 
analysis. 

FIG. 3. Phase diagram showing the pure melting line (thin line), the first-order melting with disorder (thick line) and the 
boundaries of SV pinning (dashed, Bgi, = 1.5Be) and SV creep (dotted, J = Jsb(0-8Be)); ^ = 0.8 (see text). Inset shows the 
renormalized shear modulus in units of Be, for increasing temperatures : from top to bottom, ^ = 0.45, 0.52 (critical), 1. 
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